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Characteristics of Internal Solitary Waves ' 

Garbis H. Keulegan 

This is an application of the method of approximations initiated by Boussinesq to the 
disturbances of the interface ])oints for waves of permanent form and the internal solitary 
wave. The system considered is a layer of liquid on another layer of greater density, the 
liquids of the layers being initially at rest and of constant total depth. The form of the 
wave is established. The dependence of wave velocity on wave height, on density differences 
and on layer thickness is determined. 



1. Introduction 

In investigations on the model laws of density 
currents, a project now being carried out at the 
National Hych^aulics Laboratory, need arose for the 
consideration of the genesis and the damping of 
internal waves at the interface between layers of 
fresh and saline waters. Among the many possible 
modes of such disturbances, one may consider for 
the sake of simplicity the behavior of a single 
intumescence and of nearly sinusoidal progressive 
waves. 

It has been our purpose to deal with these two 
types of internal wave motions experimentally. 
The experimental studies of sohtary waves have 
been completed. As a natural guide in the study 
of the data, resort has been made to a theoretical 
analysis. The present paper gives the basis and 
the result of this anal^^sis. Relations are here 
obtained giving the dependence of the velocity of 
wave propagation on wave heights, the form of the 
solitary wave, and the expressions for the velocity 
vector on the upper and the lower layers. Con- 
sideration of the experimental data, however, is 
reserved for a future occasion. 

After putting forth the basic conditions for the 
analysis, the question of internal waves of infinitesi- 
mal height and neghgible interfacial surface curva- 
ture is taken up in the first approximate solution. 
In the second approximate solution the character- 
istics of solitary waves are revealed. In general, 
there are seen to be some similarities between 
ordinary surface solitary waves and internal solitary 
waves. No attempt is made to extend the approxi- 
mations to a third-order analysis as no special de- 
mand is made by experimental evidence for this. 

2. Mathematical Formulation 

A layer of lighter liquid of thickness H^ and of 
density p' rests on a layer of a denser liquid of thick- 
ness H and of density p. The upper liquid at its 
free surface is exposed to air and the lower liquid 
rests on a rigid horizontal bed. Both liquids are 
initially at rest. The displacement of the interface 
with respect to its initial undisturbed position is 
denoted by h; the displacement of the free surface 
with respect to the level of the undisturbed free 
surface is denoted by A' (see fig. 1). Taking the 
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origin, 0, of rectangular axes at the undisturbed 
interface, the axis of x is drawn horizontally. The 
axis of 2 is drawn vertically and the positive branch 
points upwards. The velocity components along 
these axes are denoted by u^ w, and u\ w' , the 
primed symbols referring to particle velocities in the 
upper liquid. 

It will be assumed that the disturbances are pro- 
duced in liquids initially at rest, so that the conse- 
quent flow is irrotational and admits the velocity 
potentials </> and </>'. The vorticity that is naturally 
present at the interface will be ignored. Since the 
type of disturbance visualized is translational, that 
is, the particle velocities in vertical planes normal 
to the direction of wave motion are nearly constant 
for each layer, it is appropi'iate to introduce the 
expressions of the velocity potentials 



for the lower layer, and 



(1) 



(2) 



for the upper layer. Here, <^o, ^o, </>oj ^^^ ^o are 
functions of x and t alone. 

Altogether there are six unknowns. A., h' ^ (t>o, 0o, 
^0, and do, which are determined on the basis of the 
kinematic and dynamic conditions at the various 
boundaries. The pressure is atmospheric at the 
free surface. At the interface the pressure is con- 
tinuous. A particle in the free surface remains in 
this surface. The particles in the upper surface of 
the lower liquid at the interface remain in this 
surface. Finally, the normal particle velocities 
vanish at the horizontal bottom. 
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The mathematical formulation of the dynamic 
boundary conditions requires that expressions be 
given for the pressures at all points in the two 
layers. Assuming that the atmoshperic pressure is 
reduced to zero, 



£=^_,._l(.^+.^)+^,H' 



for the lower liquid, and 



^=^-gz~(u''+w'') + gH' 



(3) 



(4) 



for the upper liquid. Translating now the state- 
ments of the boundary conditions above into 
mathematical relations, neglecting w'^ with respect 
to u^y and w^^ with respect to ^^'^, they are 



^_,t._"...=o. 



z=H'+h', 



z = h, 
bh bh 



(5) 
(6) 



and 



at dx ' 


z = K 


, bh ,bh^ 
bt bx ' 


z=H, 


w^Q, 


z = h. 



0) 

(S) 

(9) 

(10) 



3. Internal Waves of Infinitesimal Height 

The functions ^o and 6q may be expressed in 
terms of ^q and 0o^ respectively, after noting the 
kinematic boundary conditions at the bottom and 
at the free surface. From eq 1, keeping only the 
first three terms, 



and therefore 






d</>_d^o ^^00 



bz bx dx'^ 
The latter in view of eq 10 requires that 

da: dx^ 



(11) 



(12) 



(13) 



Hence, 



<I>=4>o-(Hz+^^ 






and therefore 



b<l> 



b'4,0 



5. = -^^+^) ax- 



Similarly, from eq 2, 



<^' = </>o+^ 



bdo z^ b^<j>o 
dx 2 dx^' 



hence 



bz dx bx^ 



(14) 



(15) 



(16) 



(17) 



Taking the kinematic boundary condition of the 
upper surface, eq 6, and neglecting the product 
term which is a small quantity, 



bx 



^-iH'+h')^=. 



bh^ 
' bt' 



(18) 



And as h' may be neglected in comparison with H' , 

»2j 



hence 



be'o^ bh' „, dV; 
bx bt^ bx'' 



'^='^-'-bi^\^^-^)^^' 



and therefore 






(19) 
(20) 

(21) 



For the points at the interface, that is z=h, one 
may write from eq 14 and 15, neglecting the terms 
involving z and its square. 



and 



bt bt 



bz bx^ 



From eq 20 and 21 



and 



b^^b^ 
bt bt 



bz 



bt 



bx' 



(22) 
(23) 

(24) 
(25) 



These imply that wave height h is small in comparison 
with H or H^ and also that the curvature of the wave 
surface is small. 

Taking next the boundary conditions, eq 5, 7, 8, 
and 9, neglecting in them terms consisting of squares 
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or of products, and using the relations from eq 22 to 
25, inclusive, 






(26) 



"'(f-^O-Kt-^O' (2^) 



H 



bx^~dt' 



rj b <t>o __ 



bh' , „, b^'. 



bx^ 



bt 



-H' 



bx' 



(28) 



(29) 



Rearranging the terms in eq 27 and then differ- 
entiating twice with respect to t^ we have 



bt^ 



bt' 



Differentiating eq 27 twice with respect to x, 
b^'o _p dVo P-P' b^h 



btbx' p' btbx^ 



bx"" 



(31) 



Multiplying eq 29 by g, introducing the value of h' 
from eq 26, and differentiating the result with 
respect to t 



bx^bt 



bx^bt 



bf 



This reduces, after eliminating <^o l>v means of eq 
30 and 31, to 



-gH 



b^ 
bx'bf 



■gH'-^ 



, p b'<f>o ,„,ip-p')b^h 



p b^ , p-p' b^h 



/ bt' 



9 



be 



Differentiating the latter with respect to t, and then 
eliminating h by means of eq 28, there results, finally, 



dt' 



c)x'dt' 



which is the differential equation of wave motion in 
a liquid system consisting of two strata of dift'erent 
densities for small disturbances. 
Putting p — p'=Ap, 



o,W2=g'WH^, 
P 



and 



the wave equation, eq 32, may be rewritten 



(33) 
(34) 



/b' , b'\/b' 2 (>^\^ n ^Qs^ 



Two types of waves are possible. In one the waves 
have the velocity of propagation itcoi; in the other, 
the velocity ±co2. The plus sign refers to waves 
moving in the direction of x negative and the minus 
sign, in the opposite direction. Furthermore, these 
waves, of infinitesimal wave height and of negligible 
surface curvature, travel without distortion of form. 
The discrete values of the velocities may be shown 
from eq 34, neglecting the higher powers of the 
density differences, to be 



cot 



and 



C02 = 



'^(' + 



._ gH^H Ap / . H^H Ap 
H+H^ p V^(H+Hr p* 



(36) 



(37) 



Of these, the second refers to disturbances of the 
interface, that is, to internal waves. The first refers 
to ordinary waves; in this expression one notices the 
eft'ect of the nonhomogeneity of fluid on the velocity 
of propagation of the ordinary waves. Nonhomo- 
geneity reduces the value of the velocity of propaga- 
tion. For a given relative density difference the 
reduction is greatest when the layers are of equal 
depth. For the purposes of the present work, the 
secondary effects of the relative density differences 
will be ignored. Accordingly the velocity of propa- 
gation of internal waves of infinitesimal height, 
replacing subscript 2 by 0, will be given as 



COo = 



gH'H Ap 



(38) 



The next question to be considered is in regard to 
the particle velocities in the two layers. One com- 
mences the analysis with the lower layer. As was 
noted above, the internal waves of infinitesimal 
height progress without change of form. This fact is 
equivalent to the equality of the operators 






(39) 



The negative sign is chosen for waves moving in the 
direction of x positive and the positive sign for waves 
moving in the direction of x negative. For the sub- 
sequent analysis it will be assumed that the internal 
waves move in the direction of x positive, and thus 



dt dx 



(40) 



One now naturally selects that kinematic condition 
of the interface stating that a particle of the lower 
liquid once on the interface remains on the interface. 
This condition is given by eq 28, which now may be 
written as 



H 



^^0 



dx' 



bh 
dx' 
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and, as in terms of the potential, the particle velocity 
1^0 is 



Uo = 



(41) 



we have 



dx ~ dx 



Integrating the latter equation with respect to x and 
observing that Uq vanishes when h vanishes, then 



Uo=CO 



H 



(42) 



Accordingly, particle velocities in the lower layer are 
proportional to the velocity of wave propagation and 
vary with the wave-element height of the internal 
wave. When these heights are positive, that is when 
the wave elements are elevated, particle motion in the 
lower layer is in the direction of wave motion. 

In terms of the velocity potential the particle 
velocity in the upper layer is 



U.= ' 



bx 



(43) 



To obtain the values of Uq one now considers the 
other kinematic boundary condition of the interface. 
The reference is to eq 29, which when treated in an 
analogous manner yields 



H^Uo-j-Huo=coh[ 



(44) 



Now, when in eq 26 the partial differentiation with 
respect to t is replaced by the partial differentiation 
with respect to x, the result is 



coUo=gh\ 



(45) 



Eliminating h between eq 44 and 45, and introduc- 
ing the value of co from eq 38, one obtains the 
relation 

H'u'o(l-^jj~j+Huo=0, (46) 

and this connects the particle velocities in the two 
layers. As the term containing the relative density 
difference is small, it will be neglected, and therefore 



Uo 






(47) 



Thus the motion of the particles in the upper layer 
is oppositely directed to the motion of the particles 
in the lower layer. 

It will be instructive to show to what extent the 
internal waves do affect the free surface. This will 
be better understood if the ratio h^ jh is investigated. 
The combinations of the expressions in eq 42 and 
45 gives 



h~~ 



'gH Uo 



which, in view of eq 38 and 39, reduces to 



h~~ 



Ap H 



P H+H^ 



(48) 



The interpretation is that the displacement of the 
surface is directed oppositely to that of the interface. 
Furthermore, the disturbances of the free surface 
are very much reduced in comparison with the 
displacements of the internal waves. 

Finally, the question of the energy of internal 
waves of infinitesimal heights may be considered. 
The energy of waves is in one part potential and in 
another part kinetic. Evaluating the potential 
energy with respect to the undisturbed configuration 
of the two layers, the appropriate expression is 






Ep = g^j^ hHx + g 



2jo 



h^^dx, 



where X is the effective wave length, that is the 
length along which h^ is measm^able. Expressing 
A' in terms of h through eq 48, 

The kinetic part of the energy in terms of the 
particle velocities in the two layers is 

E,=-^^£u'o"H'dx+-^rulHdx, 

Expressing Uq in terms of iio with the use of eq 46, 
expressing Uq in terms of h and co with the use of eq 
42, and introducing the more exact values of co from 
eq 37, the final value of kinetic energy is 



['+^(ffiH-.j]j>- 






(50) 



Comparison of the expressions for the two energies 
shows that, in internal waves of infinitesimal wave 
height when the upper surface is free, the potential 
and kinetic energies are of like value. When the 
relative density-difference term is ignored the energy 
of internal waves is simply 



E=gAp hHx 



/: 



(51) 



The form of the internal waves together with the 
difference of the density of the two layers is suflBcient 
for the evaluation of the energy. 

Some of these properties of internal waves of large 
wave length and of infinitesimal wa\7e height are 
well known. A short derivation of these results, for 
example, by the method of Lord Rayleigh is given 
by Thorade [1]. 
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4. The Internal Solitary Wave 

In general when the wave-ekunent heights of the 
internal waves are finite, a deformation of the wave 
surface occurs during the travel of the waves even 
when the two liquids are assumed to be perfect and 
ideal. The extent of the deformation depends on 
the height h and on the surface curvature b^/dx^. 
The corresponding mode of the deformation in the 
ordinary translation waves, as first discussed by 
Boussinesq, is well known [2, 3]. A similar theory 
can be worked out also in the cases of internal waves, 
and this theory could form a basis for the study of 
internal waves which travel without d^ormation. A 
second procedure is to assume the existence of internal 
solitary waves, positive or negative, which travel 
long distances without deformation in the absence 
of viscosity. These intumescences will be referred 
to as internal solitary waves, and the analysis to be 
followed below will determine the form of the waves. 

Ordinarily in establishing the character of the 
progressive waves of the permanent type resort is 
made to an artifice [4]. By superposing on the flow 
a current of the magnitude equaling the velocity of 
propagation of the wave and moving in the direction 
opposite to that of the wave, the system is reduced 
to a steady state. Another method, and this will be 
used here, is that the time differentiations will be 
changed into space differentiations for the wave 
traveling in the direction of x positive using 






= — CO ^-— ; 
OX 



(40a) 



where co is the velocity of progression of the wave of 
permanent type. In the present case co is the velocity 
of progression of the solitary wave. 

The expression for the vector potential in the 
lower layer to the second approximation, is, from 
eq 1, 



</> = </>o+^ 



dx 2 dx' 6 dx^' 



(52) 



Evaluating the last term on ^i^ie right-hand side, 
since it is the term of smallest value, from the first 
approximative value. 






dx 



the vector potential now is 



dx' 



(13a) 



<^=0O+^ 



bSo 2'(>'cl>o , Hz^c)'ct>o 



dx 2 dx' 



6 dx' 



and, therefore. 






dVo , Hz^ d^c^o 



dz dx 



dx' 



2 dx^ 



In view of the bottom kinenuitic 
eq 10, 



dx c)x' 2 dx^' 



)()iiudary condition, 
(53) 



which than is the second approximative value of 
c>do/(>x. Thus, the resulting form of the vector 
potential is 



<^=<t>Q—Hz 



dVo H'z d'(t>o z' d^^o , Hz' dVo 



and, therefore, 
-H 






dx' 



dx'' 



2 dx' 2 dx'^ 



6 dx^' 



d'<t>o tP d'^o , Hz' dVo 
dx' 



2 5^4 



2 dx' 



(54) 



(55) 



This is the expression for the velocity potential 
in the lower layer to the second approximation 
involving the single unknown function </)o. When 
the height of the internal solitary wave is smaller 
than the depth H of the lower layer, the expression 
for the vertical component of the particle velocities 
at the interface may be simplified by neglecting 
z' in eq 55. Mathematically the basis of the ap- 
proximation is 



A^bVo^.^.^Vo 



2 dx 



'<^H 



dx' 



Thus, neglecting z', 
ox 






dx' 



(56) 



^^- (57) 



£)Vo 
dx' 



Now, the last two terms on t\ui riglit-liaud sitle being 
the smaUest in vahie, sul)stitutions may be made 
for them from the first ap])ioximative vahics, 



d<t>o 
dx 


-Ua= — 


d'4>o 


coo dh 


dx' 


Hdx 


d^4>o 


coo d'h 


dx^ 


Hdx' 


dVo 


COo d% 



Wo 



A 
hi 



dx* 



Hdx' 



(58) 



Introducing these in eq 55 and putting 3= A, the 
vertical component of the particle velocity at the 
interface, sign reversed, is 



dz dx' 



coS dh , o^oH' d% 



Hdx 



2 dx' 



z = h. (59) 



Similarly the expression for the vector potential 
in the upper liquid, from eq 1, is 



</> = 0o + ^ 



d^o z' d'cl>o 
dx 2 dx' 



d dx' 



(60) 
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The last term on the right-hand side being the smallest 
term, the value of ^o may be taken from the first 
approximation, 



^'L---^+H- 



bx 



bt 






(19a) 



It is permissible at this time to affect a small modi- 
fication in the analysis, and this consists of neglecting 
the term bh'/dt. The assumption implies that 

_^ ^H' ^, 
bt bx 

and may be verified as follows. In view of eq 40a, 



bh' 



bh 



~~bt-'''bi' 



Also, in view of eq 43, 45, and 38, 

bx^ dx coo c>x H Ap dx 



Thus, since 



ApHMH+H')p. 



it is permissible to neglect ?)h^/(>t and write for the 
first approximative value of bd^/dx 



bx bx^ 

Introducing this is eq 55, 

^, ^,, d^o z'b'<t>'o zm'b'<t>o 



and, therefore, 

b<i>' be'o 



bz 



da; ^ bx^ 2 ^a;*' 



(19b) 



(61) 



(62) 



Since bh' jbt may be neglected and also u'bh'/bt 
which is still smaller, the kinematic boundary condi- 
tion for the upper surface, eq 6 reduces to 



bz 
and this requires that 



z=H'+h', 



be'o b'<t>'o,H'' b% 

bx bx^ 2 da;*' 



(63) 



(64) 



which is the second approximative value of bd'o/bx. 
Substituting in eq 61, 



,._.,/ I .rr/^Vo Z' bVo,2H'' bVo zm'bVo 
0-<)oi-2/i ^^2 2 bx' 2 bx' 6 bx' 



(65) 



and 



ib£_ rj, ^_; ^^^HZ_ b^_z^ b% ,_ _. 
[d^""^ bx' ^ bx' 2 bx' 2 bx'' ^^^^ 



When the latter expression is applied to the region 
close to the interface, h' being smaller than H' , the 
term involving z' may be neglected. Thus, 

b£_ rj, b^_ ^ 5V; . HI' b^ ,. .. 

bz ^ bx' ^ bx' 2 bx'' ^^^^ 

Now the last two terms on the right-hand side being 
the smallest in value, substitution may be made for 
them from the first approximate values 



bct>o 



h 



bx~^^'~'''''H'' 



b'<l)o _ 030 bh 
bx'~H'bx 

b^ct>o^ coo b'h 
bx^ H' bx'' 



bVo 



Wo 



Hence, 



bz bx' 



_ b% 
bx' H' bx'' 



, b'<l>o coj bh cooA' b'h 



(68) 



h' bx ' 2 bx' 



z = h. (69) 



Consider next the two kinematic boundary con- 
ditions for the interface given by eq 8 and 9. These 
may be written, in view of the rule of eq 40a, 



b<h bh , bh 

bz bx bx 



z = h, (70) 



, bcb bh , ,bh , 

w'=—-^--=—o) ^ — h^ ^^^ z = h. 
bz bx bx ' 



(71) 



The last terms on the right-hand sides of these two 
equations being tha^mallest, the values of u and u' 
in them may be taken from the first approximate 
values 



U = COoYf fl^nd U'=— OJo 



A 

H' 



Thus, 



— _^— -_ blh,^ohbh 
~ bz~ ^ bx H bx 



(72) 



and 



, b(b' bh ojoh bh , ,^^, 

Substituting in eq 72 the value of w from eq 59, 
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and in eq 73 the value of w' from eq 69, the kine- 
matic boundary conditions reduce to 



H 



^Vo 



bx^ 



-0) ^-+2 -TT-^r-H — w—^r-^f 



dx 



H dx ' 2 bx^ 



(74) 



and 



-H 



, d 4>(, 



bx' 



bh woh dh ooH' b^h 

-''b^-^Wb^+-^-b^^- ^^^) 



Multiplying eq 74 by p/H, integrating with respect 
to X, and putting the constant of integration equal 
to zero since there is no disturbance at infinity, one 
obtains 



dx 



h 



P -KZ'^'^P^Jj+P^o 



a) 



\ tOpg b'h 

+ "2 bz'' 



(76) 



Multiplying eq 75 by p'jH', integrating with I'espect 
to X, and again putting the constant of integration 
equal to zero, 



, oiqH' b'h 
2 bx^' 



(77) 



"" b^^"" "H^"" ""{W) +" 
Adding eq 76 and 77, 

,b<l>'o, b(l>o r'''_L''l7>-L f-^-- ^1 ^2 



+ f[p'^'+pH]g- 



(78) 



The dynamic condition for the interface given 
by eq 7 may be written in the form 

p'f-P^=/,*-P,^+^-f- CO) 

The last two terms on the right-hand side are 
small quantities and these may be evaluated by the 
first approximative values of the particle veloci- 
ties u'=—cooh/H^ and u=oiJilH. 

Hence, 

(80) 
In view of the rule given by eq 40a, 

— p ^ ->^+ p^ >^= — yp — P )9f^ 



dx 



dx 



Expressing the value of the density difference 
p— p'=Ap in terms of ool, using eq 38, and dividing 
by CO, 



da; 






Keeping only the first two terms of the right-hand 
members of eq 54 and 65, 



and 



dx bx bx^^ 



bx bx bx"^^ 



z = hj 



z = h. 



In view of the approximative values from eq 58 
and eq 68, 



d0_d</)o , , b'^h 
bi~'b^,^'''^b^'' 

b(t>'_b(t)o, y b'^h 



z = h. 



Thus, if (cooh) b^h/bx^ is neglected, eq 82 reduces to 

"^ -bx+'bx=^''c:;iw+Hr 

+?:k^)-'(h)3 -^ «') 

Equating the right-hand members of eq 78 and 82, 
since the left-hand members are the same quantities, 
and collecting the terms, 

+1^ \w~w)^'" \m~w)\ ^' 



COo 



WH'+pH]^,=0. (84) 



If we divide by cop and ignore the terms containing 
Ap/p, the final result is 

[(~)"-][^4]-i[-i^+^]- 



1 



d^A 



This simplifies, after dividing by {H'+H)IHH', to 
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Multiplying by dh/dx^ integrating once with respect 
to X, and putting the constant of integration equal 
to zero since all disturbances vanish at infinity, we 
obtain 

Let the crest height of the wave be hi. As this is a 
maximum, c)h/()x vanishes at this point. Hence 



or 



or 



CO 

coo 



, _ H'-H , 



(87) 



V'+(w)^^- (^«) 



a; = coo 



This is the law of velocity of propagation of an inter- 
nal solitary wave. The magnitude of coo, the velocity 
of propagation when the wave is exceedingly small, 
is given in eq 38. Substituting in eq 86 from eq ^S, 



and the solution of the equation is 

h = hi sec^ ha Iy-——\ 



where 



/I {H'-H\ h, 



(89) 
(90) 



This gives the form of internal solitary waves. 

Examination of eq 86 reveals that the relation of 
the thickness of the two layers has an important 
bearing on the formation of internal solitary waves. 
When the depth of the upper layer is greater than 
the depth of the lower la3^er, H^y>H^ the internal 
solitary wave is of positive type, that is, h is positive. 
When the depth of the lower layer is greater than 
the depth of the upper layer, the internal solitary 
wave is of negative height, that is, h is negative 
everywhere. When the two layers are of the same 
depth and the difference of the densities is very small, 
the formation of a solitary wave is excluded by this 
analysis. 

The consideration of a fairly large number of 
experiments tending to verify the above theoretical 
results, first in regard to the dependence of velocity 
of propagation on wave height, and second in regard 
to the form of internal solitary waves, will be reserved 
for another occasion. 
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